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Abstract 
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1 Introduction 



Supersymmetric quantum mechanics (SQM for short) is a framework used to determine 
energy eigenvalues and eigenf unctions of quantum mechanical problems. SQM appeared 
25 years ago and has so far been considered as a new field of research, providing not only 
a supersymmetric interpretation of the Schrodinger equation, but important results to a 
variety of non-relativistic quantum mechanical problems pQ. One of the most important 
works is that Gendenshtein introduced a discrete parametrization invariance labelled " shape 
invariance" [2J- For the shape invariance potentials, their exact energy levels can be found 
analytically by making use of the supersymmetry shape invariance approach [SI IH E] • 

The association of the factorization and the hierarchy of hamiltonian method with SQM 
formalism has been introduced to obtain the approximate energy spectra of non-exactly 
solvable potentials, [01 U\ as well as the partially solvable, jBJE], the isospectral, [TP] , the 
periodic potentials, [TTj and the exponential-type potentials [I2l HSl EU HE]- Using the 
physical arguments, it is possible to make an ansatz for the superpotential which satisfies 
the Riccati equation by an effective potential. Thus, one can obtain a solution of the Riccati 
differential equation for the superpotential function. This new scheme can be successfully 
applied to obtain the spectra of the nuclear systems well fit by the deformed Woods-Saxon 
potentials. 

In this work, a study of such new exactly solvable effective potential through SQM is pre- 
sented by the hierarchy of hamiltonian method for I = and I 7^ cases. We first introduce 
SQM briefly, and then show the deformed Woods-Saxon potential hierarchy and its shape 
invariance property for the I = states. The deformed Woods-Saxon potential hierarchy is 
recovered with certain potential constants. The results obtained are improved by means of 
the shape invariance for the I 7^ states. In this case, the effective potential behaves in the 
same way as centrifugal part of the Woods-Saxon potential under the condition of the first 
approximation. Finally, this approximation is applied to the nuclear scattering process. 



2 Supersymmetric Quantum Mechanics 

The simplest way of generating a new exactly solvable Hamiltonian from a known one is 
just to consider an invertible bounded operator. In this case, a given Hamiltonian and its 
supersymmetric partner have identical spectra except the zero energy of ground state. In 
Supersymmetric Quantum Mechanics (SQM) for N = 2 we have two nilpotent operators, Q 
and Q + , satisfying the algebra 

{Q,Q + } = H S , {Q,Q} = {Q+,Q+} = 0, (1) 
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where Hs is the supersymmetric Hamiltonian. The fact that the supercharges Q and Q + 
commute with Hs is responsible for the dejeneracy. This algebra can be realized as 

where A 1 * 1 are bosonic operators. With this realization the supersymmetric Hamiltonian Hs 
is given by 

Hs= \0 A-V) = (1) H+)- (3) 
H are called supersymmetric partner Hamiltonians and share the same spectra, apart from 
the nondegenerate ground state, (see pQ for a review), 

El +) = (4) 

For the non-spontaneously broken supersymmetry this lowest level is of zero energy, E± = 0. 
We have 

where V±(r) are called partner potentials. The operators are defined in terms of the 
superpotential W(r), 

A ± = T^4- + W{r) (6) 
which satisfies the Riccati equation 

W 2 ±^=W = V ± {r) (7) 

as a consequence of the factorization of the Hamiltonians H ± . 

By definition, two partner potentials are called shape invariant if they have the same 
functional form, differing only by change of parameters, including an additive constant. In 
this case the partner potentials satisfy 

V+(r,a 1 ) = V-{r,a 2 ) + R{a 2 ), (8) 
where ai and a 2 denote a set of parameters, with a 2 being a function of a\, 

a 2 = /(en) (9) 

and R(a 2 ) is independent of r. 

Through the super-algebra, for a given Hamiltonian Hi, factorized in terms of the bosonic 
operators, it is possible to construct its hierarchy of Hamiltonians. For the general sponta- 
neously broken supersymmetric case we have 

* = -££ + *W = ^- + tf» do) 



where is the lowest eigenvalue. 

The bosonic operators are defined by (jUJ) whereas the superpotential W\(r) satisfies the 
Riccati equation 

WZ--%=W[ = V x {r)-E$\ (11) 
V 2m 

The eigenfunction for the lowest state is related to the superpotential W\ by 

^W(r) = Nexp (-^p f Q W 1 {f)df S j . (12) 
The supersymmetric partner Hamiltonian is given by 

H 2 = A^Ai + Eg* = -^-^ + [wl + -jLwi] + E£\ (13) 



2m dr 2 \ y/2m 
Thus, factorizing H 2 in terms of a new pair of bosonic operators, A 2 we get 

- A ^ + E ° 2> - + - + (14) 

(2) 

where Eq is the lowest eigenvalue of H 2 and W2 satisfy the Riccati equation, 

Wl - -4=W' 2 = V 2 (r) - 4 2) . (15) 
V 2m 

Thus a whole hierarchy of Hamiltonians can be constructed, with simple relations con- 
necting the eigenvalues and eigenfunctions of the n-members 

H n = A + n A- + Et\ (16) 

A± ^^Tr +W ^ (17) 
*?> = , ^ = (18) 

where \I/q (r) is given by (J 12)) . 



3 Deformed Woods-Saxon Potential 

We consider the following potential which is a generalization of the deformed Woods-Saxon 
potential 



l + ? e 1 a ; (l + qe 
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Here r denotes the center-of-mass distance between the projectile nucleus and the target 
nucleus. Other parameters in the potential, q is the deformation parameter (q > 1), R = 
t-qA 1 ^ 3 is the radius of the corresponding spherical nucleus or the width of the potential, A 
is the target mass number, r is the radius parameter, Vo is the the potential depth, a is the 
diffuseness of the nuclear surface and lastly C is the setting parameter which is proposed by 
us. Therefore, we can construct the hierarchy of Hamiltonians for the original Schrodinger 
equation 

h 2 

V 2 ^(r) +V(r)#(r) = £7tf(r), (20) 

2m 

and find a solution of Eq. (J2U|) by separating variables in spherical coordinates, putting 

tt(r) = - X (r)y(0,0). (21) 
r 

Then, we get the radial Schrodinger equation for all of the angular momentum states 

h 2 d 2 , . / Tl . i(i + i)h 2 \ . . 
3.1 Solution for the / = case 

For a generalization of the deformed Woods-Saxon potential given by Eq. (jl9j) . it is substi- 
tuted into the Schrodinger equation for the zero angular momentum states, 

- ^^X(r) + V(r) X (r) = E X {r). (23) 
The ground state eigenfunction Xo( r ) can be written as 

/ 2m 



Xo(r) = Nexp / Wi(r)drj , (24) 

where iV is the normalized constant. Substituting Eq. (j24|) into Eq. (}2*3~|) . we obtain 

ft iw. 

V2m, ar 

where Eq is the lowest energy-eigenvalue or the ground state energy. Through the superal- 
gebra we take superpotential 

% ( p -a(r-R ) \ 

M/ ^-7s( S ' + S ' 2 l + ,e-^» j- < 26 » 

satisfies the associated Riccati equation (Eq. lfTTjl ) and substituting this expression into 
Eq. (j23|) . we find the following identity 

h 2 S 2 ^(2^ -qg 2 ) h 2 (S 2 + aqS 2 ) (1) 



2m 2m(g + e a ( r -- R °)) 2m (o + e a ( r - fl «' 



With the comparison of the each side of the Eq. (|27|) . we obtain 

a = l/a, 



2m 
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2m 



(2S 1 S 2 - aS 2 ) = -V , 



— (S 2 2 +aqS 2 )=C. 
2m 



The eigenfunction x( r ) f° r the ground state can be expressed as 
Xo(r) = Nexp 



Ne Sir 



1 + qe- a ( r - R °1 

( p a(r-R ) \ S 2/ a 1 



e a(r-R ) _j_ g 



(29) 



Solving Eq.(|2BI> yields 



_ 2m ( C\ 1 S 2 

5l -^r Fo + 7i 2s 2 ~ v 



aq faq\ 2 2mC 

S2 = -Y ± nY) + yt- 



(30) 



At this point, using Eqs. (j2Tj) and (J30J), the supersymmetric partner potentials can be ex- 
pressed as 



v+(r) 



2m 



K_(r) = — 
2m 



S 2 l + 



g _|_ g«(r-_R ) 



+ 



si 



+ 



CX So 



aqS 2 



2Sl 



(g _|_ e a(r-R )^ q + e a(r-R ) ( g + e a(r ~ R ) 

(31) 



si 



aSo 



q _)_ e «(r-Ko) 



+ 



(g _|_ e a(r-fl )) 2 g + e^ 1 "-^) (g + e a(r--R )) 



(32) 



Clearly, one can write 



V + (r,S 2 ) = V4r,S 2 -aq) 



+ 



2m 



2m 



s. 



1 2 



9 7 25 2 2q 



2m 



2m 



w[- y ° + - 



S 2 - aq 



(33) 



q J 2(S 2 - aq) 2q 

which is precisely the requirement for the shape invariance. The shape invariant concept 
was introduced by Gendenshtein 0. In fact, comparing Eqs. (|3^|) and (JSJ) yield 



£2 —> cli , S 2 — ag — ► a 2 , 
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where the remainder R(a2) is independent of r. On repeatedly using the shape invariance 
condition Eq.(jHJ, it is then clear that 



* 2 72 h 
H(k) = + y_ {r . ^ + £ R ^ 

s=l 



2m dr 2 



(35) 



where H^ k > is a series of Hamiltonians, k=l, 2, 3,..., and if ^ = H_. Thus the bound-state 
energy spectrum of is obtained 

k 



(36) 



s=l 



and its n-th level is coincident with the ground state of the Hamiltonian H n (n = 0, 1, 2, ...). 
The energy eigenvalues of Hamiltonian are given by 



= 0, (37) 



E^ 



2m 



2m (C 



v n 



So 



2S 2 2q 



2m 



2m (C 



Vn 



2(S*2 — naq) 



S2 — naq 

2q 



(38) 

Hence, the energy levels of the deformed Woods-Saxon potential in Eq. (jl9j) for the zero 
angular momentum states are found as 

1 2 



2 r 



E n = E^ + E 



h 

2m 



2m C 



1 



h 2 { q V °) 2(S 2 - naq) 



S 2 — naq 
27 



(39) 



where note that one can write the relationship between V(r) and V-(r) is V(r) = V-.(r) + Eo- 
Substituting Eq.()3(J|l into Eq. (j3T?|) and setting C = 0, we can immediately obtain the energy 
eigenvalues 



n 2 



2ma 2 



2ma 2 V 
h 2 q(n + l) 



■ n + 1 \ 2 AmaVo 2 



h 2 q 2 



(40) 



where the minus sign of S 2 is used corresponding to q > 0. This result is in good agreement 
with the ones obtained before [TBI ITTj. 

3.2 Solution for the / 7^ case 

The Hamiltonian for the deformed Woods-Saxon potential for the I 7^ case is written as 



H = — 



2m dr 2 



1 + qe 



h 2 i{i + 1) 

2mr 2 



(41) 



The second term is taken from the deformed Woods-Saxon potential and third term comes 
from the potential barrier. The last term prevents us to build the superfamily as in the I ^ 
case, since the full potential is not exactly solvable. However, several numerical approaches 
have been utilized in order to evaluate the spectra of energy eigenvalues and eigenfunctions. 
Now, we introduce a new effective potential whose functional form is given as follows 



Vn 



h 2 l(l + l)a 2 



eff 



q + e 



2m I q + 



- R o- 



2 ■ 



(42) 



In the case of a = (1 + q)/R and for small a, the second term of Eq.(f4*2~j) behaves as a 
potential barrier term of Eq.(|41)l in first approximation and it has the advantage that the 
Schrodinger equation for this potential is solvable analytically. As the effective potential 
given by Eq. (J42|) has the same functional form as Eq.(|19|). we can solve the Schrodinger 
problem by the factorisation method of SQM and find the whole super family. Comparing 
the Eqs. (j42j) and (|19j) . one can see the transformation of C — > h 2 l(l + l)a 2 /2m. Substituting 
this parameter into Eq. (}3T)|) . we arrive at 



S 2 




+ l(l + l)a 2 . 



2 V V 2 J 

The lowest energy levels of the potential in Eq. (j41j) are given by 



(43) 



E„ 



ni 



2m 



2mVg _ l(l+i)a 
aqh 2 q 2 



+ 



l + 2n^ 

Taking a = 1/a, we can get the result as 



1 , 41(1+1) 4 

q' 2 



l + 2n + Wl + 



1 m + 1) 



(44) 



E 



h 2 



Til 



2ma? 



( 



2mVoa 2 _ 1(1+1) 
qh 2 q 2 



l + 2n+Jl + 



41(1+1) 
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1 f2mV a 2 1(1 + 1) 

2 { qh 2 lf~ 



f - | l + 2n + Jl + 



4/(/ + l) 



(45) 

This result is exactly the same with the ones obtained before for the / = states (TBI El ■ 
In our case, quantum numbers takes / = 1, 2, ... and n = 0, 1, 2, .... Thus, n = and / = 1 
corresponds to the state 2p ; n = 1 and / = 1 corresponds to the state 3p and so on. 



Now, we calculate the energy spectrum taking the deformed q = 1 and then a = 2/Rq 
or a = Ro/2. If we take q = 2 and then a = 3/Rq or a = Ro/3. For the scattering processes, 
it has been well accepted that the surface diffuseness parameter a is ranging between 0.8 
and 1.4 fm ^Hl ITH] . Here, we use the empirical value of ro = 0.90 fm, from which we only 
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retain the Woods-Saxon potential which yields "satisfactory" fits to the experimental data, 
for projectile particles and targets with mass number A. It is shown that the mass number is 
in the domains 6 < A < 30 for q = 1 and 19 < A < 101 for q = 2. Our numerical results are 
listed for the / = 1 and I — 3 in Table 1 and Table 2, respectively. They are also compared 
with exact numerical values. 

4 Conclusions 

We have applied the hierarchy of hamiltonian method in the context of SQM to get energy 
spectra of the deformed Woods-Saxon potential. We have used a new effective potential 
depending on the diffuseness parameter a in the calculations. It is a generalization of the 
deformed Woods-Saxon potential. We have obtained the exact analytical eigenfunction and 
eigenvalue for the I = case. In addition, we have also derived the solutions for I ^ 
case, using an effective potential suggested by the / = case. We have also noticed that 
a = (1 + q)/ Rq should be taken for a nuclear scattering process. Finally, we would like to 
point out that although the SQM scheme works quite well for the deformed Woods-Saxon 
potential, extensive applications to other effective Woods-Saxon-like potentials are needed 
to test the credibility of the method. 
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Table 1: Energy eigenvalues as a function of the diffuseness parameter for the state 2p 
(n = 0, 1 = 1) and q = 2. 



State 



a (fm) 



^(analytical) £ 



(numerical) 
01 



(MeV) 



2p 



0.814 
1.026 
1.148 
1.293 
1.392 



20.0 
40.0 
56.0 
80.0 
100.0 



-17.5727121 
-16.6731389 
-17.3509649 
-19.1773309 
-20.2905543 



-17.5700005 
-16.6701250 
-17.3500205 
-19.1800050 
-20.3004000 
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Table 2: Energy eigenvalues as a function of the diffuseness parameter for the states 
= 2, I = 1), 4d (n = 1, 1 = 2), 4/ (n = 0, Z = 3) and g = 1. 



States 



E 



(analytical) 
nl 



(MeV) 



E 



(numerical) 
nl 



(MeV) 



4p 



0.818 
0.969 
1.084 
1.179 



6.0 
10.0 
14.0 
18.0 



-113.754997 
-87.027172 
-73.982269 
-66.165173 



-113.7550000 
-87.0270405 
-73.9810000 
-66.1555005 



4d 



1.001 
1.109 
1.201 
1.279 



11.0 
15.0 
19.0 
23.0 



-46.633933 
-41.128693 
-37.731683 
-35.550305 



-46.6290000 
-41.1250005 
-37.7312450 



4f 



1.084 
1.201 
1.298 
1.398 



14.0 
19.0 
24.0 
30.0 



-11.191323 
-10.866104 
-10.828699 
-10.977817 



-11.1900050 
-10.8625400 
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